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Jacobi zeta function and action-angle coordinates for the pendulum 
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The action-angle coordinates for the planar pendulum problem are expressed in terms of the Jacobi 
elliptic functions and integrals. In particular, we show that the Jacobi zeta function generates the 
canonical transformation from the pendulum coordinates i? and p = d-&/dt to the action-angle 
coordinates (J, £) for both the librating pendulum and the rotating pendulum. 



I. INTRODUCTION 

It is often said that, in order to understand the Jacobi 
elliptic functions, one only needs to study the pendulum 
problem [l[ . Jacobi elliptic functions are doubly-periodic 
functions that are useful in solving many problems in 
classical mechanics, general relativity and cosmology, as 
well as the soliton solutions of several nonlinear partial 
differential equations 041] ■ For an excellent introduction 
to the Jacobi elliptic functions see Ref. [2j while a sur- 
vey of the properties of elliptic functions and integrals 
is given in Refs. @, Q- The Jacobi functions are famil- 
iar to many physicists because they transform into the 
well-known singly-periodic trigonometric (or hyperbolic 
trigonometric) functions in the limit where one of the two 
Jacobi periods becomes infinite. 

In the present paper, we show that another Jacobi el- 
liptic function, the Jacobi zeta function plays a fun- 
damental role in generating the canonical transformation 
from the pendulum coordinates •& and p = <9$ / dt to the 
action- angle coordinates (J, Q for both the librating and 
rotating motions of the pendulum. 

We begin with a brief review of the solutions of li- 
brating and rotating motion of the pendulum problem 
expressed in terms of the Jacobi elliptic functions 
The solution for the problem of the libration motion of 
a pendulum is expressed in terms of the Jacobi elliptic 
function sn(f|re) [9fl as 



$t(t, k) = 2 arcsin[-\/K sn(t|/c)], 



(1) 



where t denotes the dimensionless time and the Jacobi 
parameter k is used to define the dimensionless energy 
of the librating pendulum e = 2 k < 2. In the low- 
energy (small-amplitude) limit s < 1, we recover the 
simple-pendulum solution $i(t, k) ~ 2 yfn sini from the 
librating solution (JlJ, where sn(t|«) ~ sini. 

For the rotation motion of the pendulum (e > 2), the 
Jacobi elliptic function sn(i|«;) in Eq. (JIJ can be evaluated 
according to the relation (n > 1) 

y/Ksn(t\K) = sn (-y/K^ft -1 ) , (2) 

and Eq. ([1]) is replaced with 

i9 r (t, k) = 2 arcsin [sn (VKtjK" 1 )] . (3) 



Note that the librating pendulum oscillates between the 
turning points ±i?£o(^), where 



< i}eo( K ) = 2 arcsin(v / K) < tt, 



(4) 



while the range of motion for the rotating pendulum is 
— 7T < i? r < 7T (where — tt and tt are now considered to 
be identical points). 

The canonical momenta p = dd/dt for the libration 
solution ([TJ and the rotation solution ([3]) are expressed 



pt(t,n) = 2 \fn cn(rjft), 
p r (t,n) = 2 \[k dn {\[Kt\K _1 ) 

where we used the relations (k > 1) 

cn(t\n) = dn (^/nt]^ 1 ) 

dn(i|K) = cn (-^/k^k^ 1 ) 



(5) 
(6) 



(7) 



Using these solutions, and the trigonometric identity 1 
cosi? = 2 sin 2 (^/2), we note that the energy equation 
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(1-cosi?) = 2 k 



follows from the Jacobi identities 

cn 2 (t\K) + sn 2 (i|K) 

dn 2 (^fk~t\n- 1 



K 1 sn 2 (y/H^K x ) 



(8) 



(9) 



for the librating and rotating solutions, respectively. 

Lastly, the separatrix solution (which separates the li- 
brating solution from the rotating solution) is obtained 
by substituting the value K = 1 in Eqs. (UJ and JS|) [or 
Eqs. © and @]: 



i9 s (r.) = 2 arcsin(tanh t) 
p s (t) — 2secht 



(10) 



where we used the identities sn(i|l) = tanh t and 
cn(f|l) = dn(£|l) = scchi. Note that the turning points 
± tt of the separatrix solution (flT?]) are now reached only 
after an infinite period of time. 

The remainder of this paper is organized as follows. In 
Sec. |TTJ we introduce the action-angle coordinates (J, () 
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for each type of pendulum motion, and the librating and 
rotating solutions are reformulated in terms of Jacobi el- 
liptic functions of (J, £). In Sec. IIII| we show that the 
transformation (p, i?) — > (J, £) is a canonical transforma- 
tion and we find the generating function S for this trans- 
formation for each type of motion in Sec. IIVI Here, we 
show that the Jacobi zeta function appears naturally as 
the generating function for the canonical transformation 
(p, $) — ► (J, £) for both types of motion. 



II. ACTION-ANGLE COORDINATES 

In this Section, we calculate explicit expressions for the 
action-angle coordinates associated with the pendulum 
problem discussed in Sec. |U We begin with the action 
coordinate 



J 



— d> pdd 
2vr 1 ' 



(11) 



for the librating and rotating motions of the pendulum, 
where the magnitude of the momentum 



|p|(0,K) 



k - sin 2 (??/2) 



(12) 



vanishes at i3 = ±i9fo( K ) for the librating case (k < 1) 
while the motion of a co-rotating pendulum (p > 0) is 
separated from that of a counter-rotating pendulum {p < 
0). 

First, the action coordinate for the librating pendulum 
is calculated from Eq. (fTTj) as 



Je( K ) 



4k 



tt/2 



n J-tt/2 yl — k sir 



E(k) - (I-k)K(k) 



(13) 



where we used the substitution sin(i?/2) = y^sin^ in 
Eq. (TT2"1) , while K(k) and E(k) denote the complete elliptic 
integrals of the first and second kinds Q. In the low- 
energy limit (k <C 1), Eq. (flU)) yields J^(k) ~ 2 k. 

Next, the action coordinate for the co-rotating pendu- 
lum is calculated from Eq. (fTTj) as 



J r {n) 



2^ r' 2 

J-tt/2 



1 — k 1 sin 



(14) 



where we used the substitution — 2 <f> in Eq. (TT2")) ; the 
action of a counter-rotating pendulum (p < 0) has the 
value — J r (n) at the same energy e = 2 k. Note that 
the action coordinate is discontinuous at the separatrix 
(k = 1): J^(l) = 8/-7T and J r (l) = 4/tt. The extra fac- 
tor of 2 comes from the fact that the full librating cycle 
requires a passage from — $i$(k) to i?^o(k) and back to 
— z?£o( K )? while the full rotating cycle only requires a one- 
way passage from — 7r to 7r (which is identical to — tt). 



Once the action coordinates J(e) have been calculated 
as functions of energy e = 2 k, we can calculate the 
pendulum frequencies w(e) = {dJ/de)~ l . The librating- 
pcndulum frequency is defined from Eq. (|13p as 



, v _ (dJ t 

where we used the expression 

E(k) - K(«) 



2K(«) 



(15) 



d 

dn 



K(«). 



In the low-energy limit (k <C 1), Eq. (fT3)l yields W£ ~ 1 
(i.e., ~ 27r). The rotating-pendulum frequency, on the 
other hand, is defined from Eq. (fl4l) as 



,'dJ r 

W -( K ) = ( -37 



where we used the expression 
d r 



V^E( K - 1 ) 



K(«-i) 

K(Q 
2^ ' 



(16) 



By comparing Eqs. ([T5 ]) -([H f and Eqs. ([T3 ]) -(fTB ]) . we ob- 
tain the relations 

K^- 1 ) = VkK(k), (17) 
v / ^E(k" 1 ) = E(k) - (I-k)K(k). (18) 

Hence, the pendulum problem naturally gives the exten- 
sion of the complete elliptic integrals E and K for n > 1. 
As discussed with the separatrix solution (fit))) , the fre- 
quencies (fT5j) - (fT6|) go to zero (i.e., the periods become 
infinite) in the limit k — >■ 1 (since K goes to infinity in 
that limit). 

Lastly, the librating solutions (H} and can be re- 
expressed in terms of the action- angle (Je,Ce) coordinates 
as 



Pe{Je,(e) 



2 arcsin 



Vk sn | — Q 

7T 



- /2K 

2Vk cn — Q 

7T 



where the libration angle Q is defined as 
= u t t = (7r/2)t/K(/s), 



(19) 
(20) 

(21) 



and k(J^) is obtained from Eq. (TT3l) . In the low-energy 
limit (k < 1), Eqs. (fl9" ]> -(f20 l) become 



Pe(Je,Ce) 




(22) 



where ~ 2 /c and Q c± t. The rotating solutions ([3]) 
and ©, on the other hand, can also be re-expressed in 
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terms of the action-angle (J r ,(r) coordinates as 

tf r (J r ,Cr) = 2 arcsin sn ^K(k _1 ) ^ kT 1 ^ ,(23) 
p r {J r ,(r) = 2^dn^K( K - 1 )^ ( 24 ) 

where the rotation angle Cr is defined as 

( r = LU r t = TTy^t/KiK- 1 ), (25) 

and re(J r ) is obtained from Eq. (TT4")) . 

III. CANONICAL TRANSFORMATION 

We now show that the phase-space transformation 

(0,p) -> (C,J) (26) 

is canonical for both the librating solution [Eqs. (|19p- 
(12TJ1) ] and the rotating solution [Eqs. P3" |) -(|23 )l ] by proving 
the canonical relation 



d-& 9p 5i9 dp 



= 1 



(27) 



for each solution. 



A. Libration Case 

We first consider the case of the librating pendulum 
(re < 1). For this case (where t = C/^t)) the partial 
derivatives in Eq. (|27|) are 



d 




1 


d 




j 


LU f 


dt' 


a 




LU f 


(- 


dJ 


c 


2 


\8k 



fit dt 



(28) 
(29) 



where d/dt and djdn, are understood to be at constant re 
and t, respectively, and we have introduced the definition 



1 

^7 



m 
dn 



t 



(l-re) 



c 2re(l-re) 
We thus use Eqs. ([I9" |) -PD ]) and (J25J) to obtain 
dd 1 dd _ 2V^ 



d( out dt 
dp 1 dp 
d( LJ£ dt 



cn, 



cue 
Vre 



sn dn, 



and then, using Eq. (|29|) . we obtain 



cm 

dJ 
dp 
dJ 



2^ 



sd 

cn + 2re 



2re / dsn 1 <9sn 

dn \ dn fl# dt 

den 1 <9cn 

9re dt 



(30) 

(31) 
(32) 

(33) 
(34) 




FIG. 1: Plot of Z(t\n) versus < t < K(re) for re = 0.9999 
(solid) and Z(i 1 1) = tanht (dashed). 



where sd = sn/dn (following standard notation for Jacobi 
elliptic functions). By inserting these expressions into 
Eq. ([27]), we obtain 



d-d dp 
9C dJ 



d-d dp 
dJ dC 



a_ 

d re 



re( 



cn + sn 



= 1, (35) 



which proves the canonical relation (|2T|) for the librating 
case when the identity (|9|) is used. 

In writing Eqs. (|33 |) -(|34 |) explicitly, the partial deriva- 
tives 



<9sn 
9re 

<9cn 
<9re 



1 <9sn 
~dt 
1 den 
fTt ~dt 



Z endn 



2 re (1 - re) 
Z sn dn — re cn sn 2 



2 re (l-re) 

are expressed in terms of the Jacobi zeta function Q 



(36) 
(37) 



Z(i|re) 



dn 2 (s|re) 



ds . 



(38) 



The Jacobi zeta function (|38l) has odd parity: Z(— i|re) = 
-Z(t|re),ithasaperiodof2K(re): Z(i+2K|re) = Z(t\n), 
and it vanishes at t = nK, with n = 0,±1,±2, 
Z(nK|re) = 0. The Jacobi zeta function (I3"51) is also eval- 
uated for the separatrix case (re = 1) as 



Z(t\l) 



sech s ds = tanh t. 



(39) 



Figure [U shows the plot of Z(i|re) for re = 0.9999 (solid) 
and the plot of Z(i|l) = tanhf (dashed). In the low- 
energy limit (re <C 1), we use dn 2 (s|re) — E/K ~ i re cos 2s 
in Eq. (|3"8"]l. which yields Z(t\n) ~ \ re sin It. 

Lastly, using Eqs. ([I|-([l7]), Eqs. (jMl)-© become 



m_ 

dJ 
dp 
dl 



LU t 



2-v/k(1 - k) 
2s/k(1 - k) 



sn dn — Z ( 



(dn 2 - 



Z sn dn 



(40) 
,(41) 



and we verify again that the canonical relation (|27p is 
satisfied. 
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B. Rotation Case 



IV. GENERATING FUNCTION 



Next, we consider the case of the rotating pendulum 
(k > 1). For this case (where u = y/nt = y/n C/ W r)> the 
partial derivatives in Eq. (|27[) are 



d 




yfk~ 8 


dc 


j 


ui r du ' 


d 




uj r ( d 


dJ 


c 


T \8k 



J_ d_ 

fl r du 



(42) 
(43) 



where d/du and d/dn are understood to be at constant k 
and u, respectively, and we have introduced the definition 



f 



du 
Ok 



1 



(44) 



2(re-l) [Hk- 1 ) 

to obtain 
dn, (45) 
(46) 



We thus use Eqs. ([2"3" ]) -([2"3 | and 

80 _ 2y/K 

dp 

dc z 

where we use the overbar notation pq = pq(u|K _1 ), and 
then, using Eqs. ([2"3" ]) -([24" J) . we obtain 



— sncn, 



dd 


cu r / <9sn 


1 


<9sn\ 


dJ 


cTi \ dn 




du ) 


dp 




dn - 


h 2k 


/ddn" 


dJ 


2V« 


V dn 



_L ddn \ 

f2 r du J 



(47) 
(48) 



By inserting these expressions into Eq. (|27|). we obtain 



dfl dp dd dp d 
9C dJ ~ dJ d( = dre" 



dn 



_l 2 

re sn 



1, (49) 



which proves the canonical relation ([2"T]) for the rotating 
case when the identity ([9]) is used. 

In writing Eqs. (j4"7 ]) - ([4"8]) . the partial derivatives 



<9sn 
dn 
ddn 
dn 



1 <9sn 

fl r du 

1 9dn 

fl r du 



Z endn 



2 (re - 1) 



dnsn 



Z sncn 



2 re (re - 1 

are expressed in terms of the Jacobi zeta function [2] 
Z = ZHkT 1 ) = 



(50) 
(51) 



dn^slwT 1 ) - = ) ds, (52) 



where E = E(k _1 ) and K = H^re" 1 ). Using Eqs. ([50 
([5l|), Eqs. d47j)-((48j) become 



d>d 
dJ 
dp 
1)1 



2 (re - 1) 



Z dn — k sn cn 



(53) 



dn (ft-cn 2 ) - Zsncn ,(54) 



2y/n(K - 1) 

which again satisfies the canonical relation ([2"T]) . 



Now that we have established the canonical nature of 
the phase-space transformation ([26]) . we seek its gener- 
ating function. First, we note that the canonical relation 
([27]) can also be expressed as the differential two-form 
identity 



dp A dtf = dJ A dC, 



(55) 



where d denotes an exterior derivative (with d 2 = 0). 
The differential canonical relation (1551) therefore allows 
us to write the one-form pd"d in action-angle space as 



pdti = JdC + dS, 



(56) 



where the function J) generates the canonical trans- 
formation (|26l) . and the identity d 2 5 = guarantees the 
canonical relation (|55[) . Specifically, the generating func- 
tion S(£, J) must satisfy the partial derivatives 



dS d$ 
dC = p d( J > 
as m 
dJ = p dJ 1 



(57) 
(58) 



for the libration and rotation cases, which we now inves- 
tigate separately. 



A. Libration Case 

First, we consider the librating-pendulum case, repre- 
sented by Eqs. [p ]) -([20 ]) . ([31]) -([32]), and gOjl-gT]). Equa- 
tion (|57p for the libration case can be explicitly written 

as 



dS t 1 dS e 4k , 



d( u>£ dt cue 



cn 2 - J t , (59) 



where Eqs. ([2"0]) and (|3"T1) were used on the right side. 
Equation ([59]) can be integrated with respect to t to give 



Si(t,K) = 4k / cn 2 (s\K)ds - J e Q, (60) 
Jo 

where we assumed that Si = at t = 0. 

Next, we use the identity cm 2 = dn 2 — (1 — re), and 
the definition (13"51) for the Jacobi zeta function, to obtain 



ft t 
I Kcn 2 (s\K)ds = Z{t\K) + [E-(l-re)K] - 
Jo K 



Z(1\k) + - J e Q, 



(61) 



which, when inserted into Eq. (Ifl2"1) , yields the final ex- 
pression for the generating function 



S e (t,K) = AZ(t\n). 



(62) 
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Hence, we see that the Jacobi zeta function Z(t\n) gener- 
ates the canonical transformation (|26l) for the librating- 
pendulum case. In the low-energy limit (k <C 1), Eq. (|62|) 
yields Si ~ K sin 2 Q, where Jg ~ 2 k and ^ ~ i, which 
satisfies the conditions (|5T )) -l|58 l) with Eq. (p2| . 

Lastly, the Jacobi zeta function (f38|) has the following 
partial derivatives 



dZ 
Ok 



dZ 

dt 
1 dZ 

fk ~dt 



dn 2 — — . 



2(1 



(63) 

dnsn - Z cn) , (64) 



where fl£ is defined in Eq. (|30|) . These partial derivatives 
can be used to show that the generating function (|6"2")l 
satisfies the partial derivatives (|57| - ([55|) . 



B. Rotation Case 

Next, we consider the rotating- pendulum case, repre- 
sented by Eqs. (J23J) - (J24|) , g5])-(|46]), and ([53]) -((54]). Equa- 
tion (1571) for the rotation case can be explicitly written 

as 



dS r JH dS r 4k —2 

-ET = TT~ = — dn - J r 

uj r au uj r 



(65) 



where Eqs. (|24|) and (|45j) were used on the right side. 
Equation (|65|) can be integrated with respect to u to give 



S r (u,K) = 4\[k I dn 2 (s\k 1 ) ds — J r Cr> (66) 
Jo 

where we assumed that S r — at u = 0. 

Next, we use the definition (|38p for the Jacobi zeta 
function to obtain 



/ dn 2 (sin- 1 ) ds = Z(u\k 



Z (u\k 



E 

= u 
K 

Jr Cr 



, (67) 



where we used Eqs. (|T4|) . (|T6[) . and (|25|) . When this ex- 
pression is inserted into Eq. (|68j) . we obtain the final 
expression for the generating function for the rotation 
case 



S r (u, k) = 4\[kZ{u\k 1 ) 



(68) 



Hence, like the libration case [Eq. (|6"2")) ]. the Jacobi 
zeta function plays a fundamental role in generating the 
canonical transformation (|2T)1) for the rotation case. 

Lastly, the Jacobi zeta function (j5"2j) has the following 
partial derivatives 



dZ 



dZ 

du 

1 dZ 

Q r du 



E 

W 



du" =. (69) 
Zen - dasn) , (70) 



2k(k- 1) 



where O r is defined in Eq. (|44[) . These partial derivatives 
can be used to show that the generating function (|68|) 
satisfies the partial derivatives (|57| -(|58 l) . 

V. SUMMARY 



The problem of the motion of a pendulum represents a 
fundamental paradigm in mathematical physics. It is well 
known that its solution is intimately connected with the 
Jacobi elliptic functions, which represents an important 
class of mathematical functions that find applications in 
physics. 

In the present paper, we have shown that yet another 
Jacobi elliptic function, the Jacobi zeta function Z(i|«;), 
appears naturally in the canonical transformation (|26[) 
that defines the action-angle coordinates for the pendu- 
lum problem. Indeed, it is used to generate the canonical 
transformation for the libration motion (|62j) and the ro- 
tation motion (1551) of the pendulum. 
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